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Abstract. We present the analysis of deep colour-magnitude diagrams (CMDs) of 6 stellar fields in the LMC. The 
data were obtained using HST/WFPC2 in the F814W (~ I) and F555W (~ V) filters, reaching V555 ~ 26.5. We 
discuss and apply a method of correcting CMDs for photometric incompleteness. A method to generate artificial 
CMDs based on a model star formation history is also developed. This method incorporates photometric error 
effects, unresolved binaries, reddening and allows use of different forms of the initial mass function and of the SFH 
itself. We use the Partial Models Method, as presented by Gallart and others, for CMD modelling, and include 
control experiments to prove its validity in a search for constraints on the Large Magellanic Cloud star formation 
history in different regions. Reliable star formation histories for each field are recovered by this method. In all 
fields, a gap in star formation with r ~ 700 Myrs is observed. Field-to-field variations have also been observed. 
The two fields near the LMC bar present some significant star forming events, having formed both young (r 
~ 1 Gyr) and old (r 10 Gyr) stars, with a clear gap from 3 — 6 Gyrs. Two other fields display quite similar 
SFHs, with increased star formation having taken place at r ~ 2 — SGyrs and 6 r ~ lOGyrs. The remaining 
two fields present star formation histories closer to uniform, with no clear event of enhanced star formation. 

Key words, galaxies: galaxies: formation, galaxies: evolution; stars: statistics 



1. Introduction 

The study of the star formation history (SFH) in near 
galaxies can contribute to our understanding of the evo- 
lution of these objects, the mechanisms involved in the 
star formation bursts, their associated time scales, and 
the initial stellar mass function (IMF). Star formation 
time scales and IMFs are particularly important, being 
used as input for population synthesis models which try 
to interpret integrated light from distant galaxies to ob- 
tain information about their component stars. Moreover, 
by matching SFH to dynamical models, it is possible to 
investigate the influence of gravitational interactions on 
star formation. 

The Large Magellanic Cloud (LMC) is a satellite of 
the Galaxy, being at the same time quite distinct from it 
and yet near enough to it to allow a study of its structure, 
internal kinematics and stellar populations (Westerlund 
1990). By its proximity and nature, the LMC provides an 
excellent laboratory for the study of the processes involved 
in star formation, such as interactions with neighboring 
galaxies (especially the Galaxy and the Small Magellanic 
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Cloud (SMC)), the influence of bar dynamics and chemical 
enrichment (Gardiner et al. 1998). 

The main advantage of LMC stellar population stud- 
ies is the ability to build detailed colour-magnitude dia- 
grams (CMDs) in various regions within it. Furthermore, 
the LMC supplies a unique chance to compare the inferred 
SFH through the study of held stars with the SFH found 
from studies of its cluster system. The reason is that the 
LMC presents a rich cluster population of largely varying 
ages. Cluster CMDs have the advantage of being com- 
posed of a single population, making their analysis much 
simpler. Studies using LMC clusters show that the major- 
ity of them are relatively young, with ages t ^ 4 Gyr, a 
small number has r 10 Gyr and only one is of interme- 
diate age (r ~ 8 Gyr) (Mateo, Hodge & Schommer 1986). 
However, only a small fraction of LMC stars belong to 
clusters, rendering the held CMD analyses of more gen- 
eral interest and more representative of the global SFH. 

Ground-based CMD studies have favoured a predom- 
inance of young or intermediate-age populations in the 
LMC. Bertelli et al. (1992), for example, have analyzed 
CMDs in a few LMC outer regions and modelled them 
with a single-burst SFH with r ~ 2 - 4 Gyr, with a factor 
of ^ 10 increase relative to the average quiescent star for- 
mation, the exact burst age being dependent on the details 
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of stellar evolutionary theories available then. The same 
group also found some evidence for field-to-field variations 
in SFH (Vallenari et al. 1996). 

Earlier ground-based work did not reach deep enough 
magnitudes to directly access the main sequence turn-offs 
of stellar populations with ages t ^ 6 Gyr. This became 
possible only with the use of the Hubble Space Telescope 
(HST), whose images were free of atmospheric turbulence 
and, therefore, much sharper and deeper. Gallagher et al. 
(1996), based on a single HST Wide Field and Planetary 
Camera 2 (WFPC2) field, suggested a roughly constant 
star formation rate (SFR) in the past few Gyr, except for 
a small period of enhanced star formation about 2 Gyr 
ago. Holtzman et al. (1997) analyzed the luminosity func- 
tion of the same HST field, located some 4° from LMC 
bar, and favoured a constant SFR for 10 Gyr in the LMC, 
followed by an increase factor of 3 in the last 2 Gyr, yield- 
ing roughly the same number of stars younger and older 
than 4 Gyr. They also identified an intermediate-age turn- 
off in their field CMD. Elson, Gilmore & Santiago (1997) 
studied an HST field closer to the LMC bar and identi- 
fied two main populations, with ages r ~ 1 - 2 Gyr and 
r ~ 2 - 4 Gyr, the exact ages depending on their metal- 
Hcity. Geha et al. (1998) have studied 3 HST fields in the 
outer region of the LMC and also suggested a model with 
roughly equal numbers of stars younger and older than 
r = 4 Gyr. Studying fields inside the LMC bar, Olsen 
(1999) suggested that the star formation has been more 
continuous, possibly extending back for a longer period 
and taking place at a roughly constant SFR as compared 
to the outer fields. This differs from the conclusions of 
Elson, Gilmore & Santiago (1997). Holtzman et al. (1999) 
have not found obvious evidence for bursts of star forma- 
tion and concluded that the field SFH differs from that 
based on the LMC cluster age distribution. 

In this paper, we describe the CMD modelling of 6 deep 
HST fields located in different directions along the LMC, 
all of them within 6° of the LMC centre. Our main goal is 
to recover information about the SFH of the LMC field po- 
pulation, the age-metallicity relation Z{t) and IMF. Our 
samples have been combined with those of a previous pa- 
per (Castro et al. 2001) and make up an extended, deep 
and homogeneous photometric dataset. The paper is out- 
fined as follows: in Sect. 2 we describe the HST/WFPC2 
data, the stellar samples derived from them, and the pho- 
tometry, including methods of accounting for sample in- 
completeness. The combined sample is presented in Sect. 
3; in Sect. 4 we present the CMD modelling algorithm, 
which is essentially the Partial Models Method (PMM) by 
Gallart et al. (1999), Aparicio et al. (2001) and Carrera 
et al. (2002), and discuss the statistical tools used in our 
model vs. data comparisons; finally, in Sect. 5 we present 
the results, which are discussed in Sect. 6. 

2. Observations and data reduction 

We use data obtained with the Wide Field and Planetary 
Camera 2 (WFPC2) on board HST for 6 fields near 



Table 1. Basic characteristics of the our HST/WFPC2 
fields. The columns show the field name, equatorial coor- 
dinates (J2000), angular distance with respect to the LMC 
centre, total exposure times and the number of exposures 
in I and V. Note that the identification of each field is 
related to the rich LMC cluster parallel to which it was 
observed. 



Field 


a (h m s) 


5 (deg ' ") 


Ang. 


ti 


tv 


Ni 






(2000) 


(2000) 


Dist. (°) 


is) 


(s) 






NGC1805 


5 01 42 


-65 59 58 


4.28 


2200 


2200 


2 


2 


NGC1818 


5 05 00 


-66 25 13 


3.75 


4800 


7200 


4 


6 


NGC1831 


5 05 21 


-64 49 50 


1.51 


2200 


2200 


2 


2 


NGC1868 


5 13 48 


-63 52 10 


5.96 


2200 


2200 


2 


2 


NGC2209 


6 07 39 


-73 46 07 


5.28 


4800 


7200 


4 


6 


Hodge 11 


6 15 07 


-69 49 08 


4.44 


4800 


7200 


4 


6 



the rich LMC clusters NGC1805, NGC1818, NGC1831, 
NGC1868, NGC2209 and Hodge 11. These data are part 
of the GO7307 project entitled "Formation and Evolution 
of Rich Star Clusters in the LMC" (Beaulicu et al. 1999; 
Kerber et al. 2002; de Grijs et al. 2002). For each field, a 
set of at least 2 exposures were obtained using each of the 
F814W (I) and F555W (V) filters. The main properties of 
the fields (hereafter identified by their nearby rich LMC 
cluster) are given in Table ^ columns 2 and 3 list their 
equatorial coordinates (J2000), followed by their angular 
distance from the LMC optical centre, total I and V expo- 
sure times and the number of individual exposures taken 
with each filter. 

All exposures were put through the standard HST 
pipeline procedure that corrects them for several instru- 
mental effects, such as bias and dark currents and ffat- 
fielding (Holtzman et al. 1995a). The exposures taken in 
each field/filter configuration were then combined using 
the IRAF'^ task CRREJ, in order to increase the signal to 
noise ratio and remove cosmic rays. The final combined 
image at each configuration was then used for sample se- 
lection and photometry. 

2.1. Sample selection 

The IRAF DAOPHOT package was used to detect and 
classify sources automatically in each combined image, i. e. 
the sample selection was carried out independently in each 
photometric band. Star candidates were detected using 
DAOFIND, with a peak intensity threshold for detection 

^ IRAF is the Image Reduction and Analysis Facility, a gen- 
eral purpose software system for the reduction and analysis 
of astronomical data. IRAF is written and supported by the 
IRAF programming group at the National Optical Astronomy 
Observatories (NOAO) in Tucson, Arizona. NOAO is operated 
by the Association of Universities for Research in Astronomy 
(AURA), Inc. under cooperative agreement with the National 
Science Foundation 
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set to 5 CT, where cr corresponds to the rms fluctuation in 
the sky counts, determined individually for each combined 
image. 
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Fig. 1. Observed colour-magnitude diagram of the field 
near NGC1818. 

A preliminary aperture {radius — 2") photometry was 
carried out by running the task PHOT on all objects out- 
put by DAOFIND. Star/galaxy separation required mod- 
elling point sources by fitting the profiles of bright, isolated 
and non-saturated stars to a Moffat function (/3 = 1.5). 
We used the IRAF task PSF for that purpose. This PSF 
model was then fitted to all objects found by DAOFIND 
with the ALLSTAR task. 

The output parameters given by ALLSTAR {x^, 
sharpness, magnitude uncertainties 61 and SY) are useful 
to separate a purely stellar sample from the remaining de- 
tected sources, like distant galaxies. Stellar sources present 
small 61 and 6Y values. Thus we threw away objects with 
61 or 6V greater than a constant value, typically ^ 0.15 
mag. These cut off values were chosen separately for each 
field/filter, based on visual inspection of the more conspic- 
uous stellar and non-stellar cases. This stellar/non-stellar 
classification was checked by visual inspection of the stel- 
lar sample selected and iterated when necessary, until a 
"clean" stellar sample resulted. The final stellar sample in 
each field, and which was used in the subsequent CMD 
analysis and modelling, was the merging result of stars 
selected in both bands, i. e., with small values of ^I and 
6Y measured by ALLSTAR. 

Figuren]shows the final observed CMD of the field near 
the cluster NGC1818. One can see a main sequence (MS) 
that stretches down to V555 ~ V = 26.5 . The saturation 
limit, also found using the ALLSTAR 61 and 6V parame- 
ters (saturated stars display large magnitude uncertainties 



from PSF fitting), was found to be about V555 = 18.5 for 
this field. The red giant branch (RGB) is visible as well, 
along with a turn-off around V555 ~ 22.5 This latter cor- 
responds to an old population (r ^ 10 Gyr, as shown in 
Castro et al. 2001). The red giant clump (RC), found in 
intermediate-age populations, is also visible at V555 ~ 19 
and V555 — I8i4 ^ 0.9 . The faint red stars at the right 
of the CMD are likely low-luminosity M dwarfs in the 
Galaxy. The other 5 deep LMC fields in our sample dis- 
play the same general features described above. 

The data suffer from an important effect: sample in- 
completeness. Our CMD modelling algorithm (presented 
in §4) therefore needs to incorporate such an effect in or- 
der to place models and data on equal footing. Quantifying 
photometric uncertainties and applying them to the model 
CMDs is extremely important as well, as they are respon- 
sible for most of the observed CMD spread. These data 
corrections are the subject of the next subsections. 

2.2. Photometric uncertainties 

There are two approaches to quantify photometric uncer- 
tainties: an empirical approach, through the comparison 
of independent magnitude measurements of the same ob- 
ject, or a semi-empirical approach, in which the photo- 
metric uncertainties follow from the measured signal and 
an adopted model for the associated noise, such as that 
used by the IRAF task PHOT. In the empirical approach, 
one needs two or more images of identical characteristics, 
such as total exposure time, detector noise, etc. We un- 
fortunately do not have two identical combined images for 
each field and filter. 

We tried to use the independent magnitude measure- 
ments taken from individual exposures of each filter/field 
configuration in order to empirically estimate the photo- 
metric uncertainties in the combined images. However, the 
photometric error measured in a single exposure tends to 
be greater than that of the combined image, due to the 
smaller signal to noise ratio. Assuming Poissonic statistics 
of photon counts, overlooking detector noise and assum- 
ing N individual exposures of the same exposure time, 
the scaled magnitude uncertainty of the combined image 
should be: 



6m I 



1 



(1) 



However, we failed to confirm the scaling in mag- 
nitude uncertainties given above, which suggests that 
the assumptions that lead to eq. ^ are over-simplistic. 
Therefore, we decided to rely on the model uncertainty 
from the IRAF PHOT task. These uncertainties take de- 
tector (readout) noise and background (sky) noise into ac- 
count as well. The spread of stars in the CMD depends on 
the photometric uncertainties of both bands. This means 
that the expected spread in colour cTcoiour is given by 



2 2 

^colour = 



(2) 
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2.3. Sample completeness 

In any image, bright stars, which are not the majority, 
are easily detected. The faint ones are often not detected. 
Photometric completeness is defined as the percentage of 
the total number of stars that were successfully included 
in the sample by the detection and selection process. The 
knowledge of the completeness function is essential when 
the objective is to compare an observed CMD to model 
ones, with the goal of reconstructing an SFH. We thus car- 
ried out experiments to estimate the completeness func- 
tion for each of our HST/WFPC2 fields. These experi- 
ments involved adding the same artificial stars, with a 
given magnitude and color, to the combined V and I im- 
ages, and submitting these images, with real plus artificial 
stars, to the same detection and selection process as the 
real data. Thus, the extra stars were subjected to auto- 
matic detection with DAOFIND, to aperture photometry 
using PHOT, to PSF fitting and star/no-star separation 
using ALLSTAR and, finally, to the final matching of stel- 
lar lists in the I and V bands, as described in Sect. 2.1. 
Most previous works estimate completeness as a single 
variable function, such as C(V) or C(I), often overlook- 
ing variations in completeness along the color axis in a 
CMD. Others have computed single variable functions in- 
dependently and have then taken their product in order to 
quantify completeness at any given CMD position (Mateo 
1988). However, since the position within the field (and 
therefore the proximity to a bright star, for instance) is 
the same in both filters, the detection of a given star in 
the combined V image is not an event independent of de- 
tection of the same star in the combined I image. Thus, 
this latter approach invariably leads to an underestimate 
of completeness, the amplitude of the bias depending on 
issues such as field crowding, exposure times, etc. Our ap- 
proach, by assigning a magnitude and colour to each ar- 
tificial star and keeping the same CCD position for both 
filters, more effectively reproduces the true bi-variate com- 
pleteness function. This resulting completeness function is 
then expressed as: 

, , recovered artificial stars in both V and I , . 

C(V 1, V) = ; — — ; — (o) 

input artificial stars 

Figure [3 shows an example of this bi-variate complete- 
ness function C(V — I,V). The completeness values are 
shown in percentage values (or C(V — I, V) x 100%). They 
are superposed on the CMD itself. It is possible to ob- 
serve that, at faint magnitude levels, C(V — I, V) decreases 
with V magnitude, as expected. For V ^ 26, however, the 
data seem to be largely complete. Moreover, C(V — I, V) 
is larger in the redder part of the low main sequence (for V 
XL 26). The latter effect may be easily understood as part 
of the selection process: faint stars in V, but bright in I 
(i.e., redder stars) have a larger probability to be selected 
than faint stars in both bands. This effect yields larger 
C(V — I, V) for redder V-I colours. 

The bi-variate C(V — I, V) function is used as a prob- 
ability that a model star, generated according to the al- 
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Fig. 2. Observed colour-magnitude diagram of the field 
near Hodge 11 (dots). The number grid shows the 
C(V — I, V) X 100% completeness values for the same 
field, calculated according to the method described in §2.3. 

gorithm described in §4.1, is actually included in the sam- 
ple. The details are left for that section. It is clear that 
the completeness function values suffer from some uncer- 
tainty. Therefore, our CMD analysis shown in §5 was re- 
stricted to CMD regions where completeness was close to 
unity (typically V ~ 24.5 — 25.0). Notice that the lower 
parts of the CMD are not crucial to the SFH reconstruc- 
tion process, since the SFH is more strongly influenced by 
main-sequence turn-offs and the positions of evolved stars. 

3. Combining different samples 

In this work we have done the sample selection and pho- 
tometry for 6 LMC fields. In Castro et al. (2001), 7 inde- 
pendent (or almost) LMC fields were studied. The data 
from both studies have a common origin, since they were 
observed with the same instrument and as part of the 
same project. All these fields are parallel HST/WFPC2 
observations of a rich LMC cluster and therefore lie about 
7.3 arcmin from this cluster. Therefore, our 6 fields can be 
paired with 6 of the fields studied by Castro et al. (2001). 
We thus decided to combine the stellar samples in each 
pair of fields in order to enhance the number of stars used 
in the CMD analysis. 

In Figure 13 we compare the stellar luminosity func- 
tion (LF) (upper panel) and the CMD (lower panel) of 
the two fields close to the rich LMC cluster NGC1831. 
The LF from our sample (dotted line in panel 3a) has 
been corrected for completeness, as described in §2.3. The 
LFs are very similar, the differences in star counts being in 
general consistent with random fluctuations (a/ZV bars are 
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shown for the Castro et ah data). At the brightest magni- 
tude bins (V555 ^ 19.5), the excess of stars from Castro et 
aL reflects their brighter saturation limit, caused by their 
shorter exposure time. At intermediate magnitudes, there 
may also exist a slight systematic excess of stars in the 
field studied by in Castro et al. (2001), suggesting that 
field-to-field variations in number counts are not entirely 
random. Even though the Castro et al. LF has not been 
corrected for sample incompleteness, this effect manifests 
itself only fainter than V555 ~ 25. The CMDs (panel 3b) 
from both samples match one another quite naturally, de- 
lineating the same loci. The RGBs are slightly displaced 
from each other (by A(V555 — /814) ^ 0.05, in the sense of 
our data (open circles) being redder). This behavior is not 
uncommon in the other fields as well and probably result 
from PSF variations with time, which may in turn reflect 
optical variations, such as changes in telescope focus. 

In order to place both samples in a common photomet- 
ric system and thus avoid artificial spread in the combined 
sample CMD, we apply colour offsets to Castro et al. data. 
This was done by first defining a fiducial line for each sep- 
arate sample, formed by the mean colour values at 0.25 
magnitude bins along the MS. We then changed the /814 
magnitude values of Castro et al. (2001) so that the two 
fiducial lines would match each other. For the RGBs, we 
found a single colour difference and changed the I814 mag- 
nitudes from Castro et al. to bring them together. Typical 
corrections in colour were of ~ 0.05. Notice that this ap- 
proach, besides being simple, yields a homogeneous pho- 
tometric dataset, while preserving as much as possible of 
the original (V555, V555 - Ism) CMD. 

The final combined sample had to necessarily respect 
the differences in exposure time, sample selection and 
completeness corrections. Exposures times of Castro et 
al. 2001 data were smaller than ours. Consequently, their 
CMDs have more information on the bright part of the 
CMD, while our data have tend to better sample the faint 
part. Also, they applied no completeness correction to 
their sample. Therefore, we defined the following combi- 
nation scheme: 

1. For the bright CMD region (V < 19.5) only the 
Castro et al. stars were used; 

2. For the faint CMD (V > 23.5) — > only our stars were 
included; 

3. For the intermediate CMD (19.5 < V < 23.5) both 
samples were included; 

The upper magnitude limit for our sample was meant 
to conservatively avoid any effect caused by saturation in 
our photometry. Likewise, the lower limit for Castro et al. 
stars is a conservative estimate of their completeness mag- 
nitude limits. The combination scheme adopted, of course, 
significantly alters the number counts along the observed 
CMD, since the effective solid angle covered by the sam- 
ple is twice the size in the intermediate CMD regions than 
in the upper or lower ones. In Sect. 4 we describe how to 
incorporate this sampling effect in a model CMD. This is 



Table 2. The result of combining the two photometric 
samples as discussed in §3. The columns show, respec- 
tively, field name, original number of stars in our sample 
and number of stars in the combined one. 



Field 


T^stars 


#stars 




our sample 


composed sample 


NGC1805 


2576 


3857 


NGC1818 


3566 


5578 


NGC1831 


1040 


1623 


NGC1868 


806 


1395 


NGC2209 


1019 


1483 


Hodge 11 


1420 


2143 



done simultaneously to the incorporation of completeness 
effects. 

In two fields, the ones close to NGC1805 and 
NGC1831, some positional overlap exists between the 
fields studied here and those by Castro et al. (2001). 
Duplication of stars was prevented by using only our data 
in these overlapping regions. Positional matches between 
the two samples were found using the astrometry from 
the task METRIC available in the IRAF.STSDAS pack- 
age. As the astrometric solutions for the overlapping fields 
were not unique, small positional offsets (~ 0.5") were ap- 
plied to the Castro et al. positions in order to make them 
comparable to ours. 

Table El shows the results of our sample combination: 
column 2 lists the number of stars in our sample alone, 
whereas column 3 lists the number of stars in the com- 
posite sample. Both numbers correspond to stars brighter 
than the adopted completeness limit V ~ 24.5 — 25.0. The 
final stellar sample is markably greater than our original 
one, which translates into better statistics and more reli- 
able CMD modelling. 

Figure^lshows the CMDs of the final combined sample 
for all 6 fields, which are again identified by the rich cluster 
close to which they lie. These may be directly compared 
to the CMDs shown by Castro et al. (2001) (figures 3 and 
4). Apart from variations in magnitude limits, our final 
CMDs display the same features and have similar MS and 
RGB positions and widths to those of the previous work. 
In particular, the RGB in the NGC1805 region remains 
much broader than the other RGBs, something that was 
interpreted as due to differential reddening. Some can- 
didate turn-offs of intermediate populations identified by 
Castro et al. (2001) still remain, while others have proven 
to be artefacts caused by limited star counts. 

4. CMD modelling and statistical tools 

4.1. Synthetic CMD 

CMDs have precious information concerning the star for- 
mation history of a given population. From observed 
CMDs one can identify MS turn-offs of various ages, there- 
fore constraining the number of star forming bursts and 
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Fig. 3. (a) Luminosity functions of the field near 
NGC1831. Our sample is shown as a dotted line and 

that of Castro et al. (2001) as a solid line with Poissonic 
bars, (b) Observed colour-magnitude diagrams of our 
NGC1831 field (open circles) and that of Castro et al. 
(2001) (crosses). 



their ages, place limits on the youngest stars from MS ter- 
mination and also constrain the metal enrichment, based 
on isochrone fitting to turn-offs and the RGB. Star counts 
in selected CMD areas are also usehil in discriminating 
star formation scenarios. A viable technique, resulting 
from ever improving stellar evolution models and data 
quality, is based on building theoretical CMDs according 
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Fig. 4. CMDs of the final combined samples for each of 
the fields, which are identified by the rich cluster close to 
which they lie. 



to a full model SFH. By comparing observed CMDs to 
synthetic model ones, we can then obtain constraints on 
the SFH that make use of all the information available. 
The process of generating artificial stars and building the 
synthetic CMD used in this work is as follows: 

1. We choose a model SFH, by assigning relative star for- 
mation rates (SFRs) as a function of time, spanning 
the range of ages to be covered, typically 8.00 ^ logir) 

< 10.25. 

2. We also adopt an age-metallicity relation, Z{t). In the 
present analysis we chose: 

- Z{t) = 0.008 for T < 950 Myrs; 

- zIt) = 0.006 for 1 Gyr ; 

- zIt) = 0.004 for 2.2 Gyr < r 

- Z{t) = 0.002 for 7.6 Gyr < t < 12.6 Gyr and 

- Z{t) = 0.001 for r > 12.6 Gyr 

This choice of Z{t) is representative of the available 
data on LMC clusters (Olszewski et al. 1996). 

3. We use an isochrone set that covers our age and 
metallicity intervals. In this paper we have used the 
isochrones from Girardi et al. (2000); Given a SFH 



< 2.1 Gyr; 
^ 7 Gyr and 
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{SFR{t)), each isochrone contributes with Ntheo,i 
stars, given by: 

NtheoA-ri-Ti) ^ C SFR{t) ^{m)dmdT{4) 

where the age hmits ti and T2 are chosen from the in- 
tervals in the isochrone grid (in the case of Girardi et 
ah 2000, successive isochrones are placed Alogr — 0.05 
apart). ip{m) is the chosen IMF and the limits in the 
mass integral are those that map onto the magnitude 
limits of the observed CMD, using the mass-luminosity 
relation embedded in the isochrones themselves. For 
the IMF we used that of Kroupa (2002), there being 
currently little variation in IMF shape inferred for dif- 
ferent populations (Kroupa 2002). The mass limits will 
depend not only on the Padova stellar models, but also 
on the adopted distance modulus, since our apparent 
magnitude CMD limits will map onto different abso- 
lute magnitude limits as this latter parameter is varied. 
The intrinsic distance modulus at each field position 
was allowed to vary around a reference value taken 
from a model for the LMC disk inclined by 45° rela- 
tive to the line of sight, with the line of nodes aligned 
in the north-south direction and (m — M)o = 18.5 for 
the LMC centre (Westerlund 1990, see also Castro et 
al. 2001 and §5 for more details). C is a constant of 
normalization, chosen so as to match the chosen total 
number of model stars (given by J2i ^theo,i{Ti,T2)). 
4. Having found the number Ntheo.i of stars to be drawn 



frt 



the 



nth 



isochrone, we randomly draw stellar 



masses according to the IMF of Kroupa (2002) and 
compute the corresponding I and V magnitudes given 
by the isochrone and the chosen intrinsic (m — M)o. 
The masses are drawn as follows: given a random num- 
ber r, between and 1, we solve the equation below 
for m: 



ip{m)dm — r 



ip{m)dmdT 



(5) 



5. We allow for unresolved binarism by repeating step (4), 
assuming that all binaries are coeval. We then combine 
the I and V fluxes of the two components and compute 
the corresponding system magnitude and colour. As 
the binary fraction, /(,i„, is largely unconstrained for 
LMC field stars, we explored different values of fbin 
and tested the sensitivity of our models to different 
choices of this parameter (see §5). 

6. We then apply the reddening vector (Av,E(B — V)) 
to the system magnitudes and colours, defining its 
theoretical CMD position. For this purpose we use 

3.1 and the photometric transformation to 



E(B-V) 

the vegamag WFPC2 system according to Holtzman 
et al. (1995a); the E(B — V) values primarily used were 
taken from Castro et al. (2001), although we also tried 
alternative values of this parameter (see §5). 



7. We spread our star positions in the CMD based on 
a Gaussian distribution of errors with try and ai as 
measured from PHOT/IRAF (see §2.2). 

8. We verify whether the star (or binary system) ends 
up inside the magnitude limits of observed CMD and 
discard it if it does not. 

9. With the observed CMD position at hand, we then 
incorporate the sampling effects. A probability p of 
actually getting into the CMD is assigned to each star 
and compared to a randomly chosen number r. If r < 
p the star is included, otherwise it is discarded. The 
probability is p = 0.5 C(V - I, V) in the CMD regions 
where only one of the photometric samples (either ours 
or from Castro et al.) was used, and p = C(V — I, V) 
in the CMD regions where both samples contributed 
to the observed CMD. 

In Figure[Slwe show some examples of synthetic CMDs: 
(a) a pure old population (r 9 Gyr); (b) a pure 
intermediate-age population (r ^ 2 - 3 Gyr); (c) a pre- 
dominantly young population (r ^ 1 Gyr) and (d) a mix- 
ture of an intermediate-age and an old population (r ~ 1 
- 3 Gyr and t ^ 6 Gyr). In figure ISJa) we see only old 
stars and, as a consequence, the CMD shows a clear single 
turn-off whose magnitude V555 is fainter and redder than 
those seen in the other panels. Moreover panel (a) shows 
a horizontal branch, common to an old population CMD. 
A red clump is visible in panels (b) and (d), as a conse- 
quence of intermediate-age stars. Panel (c) is the only one 
whose main sequence stars extends towards the brightest 
part of CMD (V555 ^ 19), with very few evolved stars. 

4.2. Statistical tools 

For CMD modelling we used a direct comparison of the 
distribution of stars in the observed CMD with model 
CMDs resulting from a large number of possible SFHs. 
The comparison method used was the Partial Models 
Method (PMM), developed by Gallart, Aparicio and 
collaborators and described in detail by Gallart et al. 
(1999). Other methods that reconstruct a SFH from direct 
CMD comparisons do exist (Hernandez, Valls-Gabaud & 
Gilmore 1999, Dolphin 2002), some of which have been 
applied to LMC fields. Our choice of PMM was based on 
the simplicity of implementing it. 

The approach in the PMM is to create several synthetic 
CMDs from a constant SFR{t), each CMD covering an 
age interval; these CMDs are called partial models. CMDs 
from any function SFR(t) can then be defined as a linear 
combination of the partial models. If we divide the CMD 
into a given number of regions or boxes, the number of 
model stars expected to be found in the j*^ CMD region 
is given by: 



i=l 



(6) 



where m is the number of partial models used, iV^.j 
is the number of stars from partial model i located in the 
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(b) 



> 22 
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"555 ^814 



Fig. 5. Some synthetic CMDs. (a) a pure old popula- 
tion; (b) a pure intermediate-age population; (c) a pre- 
dominantly young population and (d) a mixture of an 
intermediate-age and an old population. 

CMD region mentioned. The set coefficients can be var- 
ied arbitrarily, accounting for any shape for the SFR{t), 
with the only constraint that 



= 1 



(7) 



The coefficients are relative SFR values, i.e., a/ < 
a >= SFR/ < SFR >. The partial models are created 
from a constant SFH and by fixing the normalization con- 
stant C in eq. 0]for all age intervals [ri, T2]. To obtain the 
SFH that best describes the data CMD, we use essentially 
the xt statistical tool defined by Gallart et al. (1999), 
but with the modification proposed by Mighell (1999) and 
Dolphin (2002), which is shown below: 



and 



1^ [ffik^Ml^ 



(8) 



= r - 1 (9) 

In eq. |H1 Nqj is the number of observed stars in the 
j*'* CMD region , r is the number of CMD regions and / 
is the normalization factor, which is the ratio of the to- 
tal numbers of observed to model stars. The xt parame- 
ter allows one to discriminate, among the different SFHs 
(represented by different sets of a^), those that best re- 
produce the observed data. We have to keep in mind the 
non-uniqueness of the best solutions. One needs a criterion 



to select satisfactory models. We defined the acceptable 
range of xt values as being 



ixDrnin < xl < ixDrnin + 2(75 



(10) 



where as represents the standard deviation from {xt)min 
of 100 realizations comparing the best model CMD (the 
solution that yields {xt)min) to the observed CMD. We 
considered the average values among these solutions as 
representing the best SFH, using the scatter around this 
average as an uncertainty estimate. 

SET A 




Fig. 6. Set A - artificial CMD composed by 8 partial mod- 
els. Age intervals are showed in Gyrs. Each panel shows 
two partial models. The regions used in the modelling pro- 
cess are indicated. 

The choice of the number of CMD regions in the data 
X model comparison was driven by the need to sample 
different CMD positions while keeping a significant num- 
ber of stars in each region. We defined r = 13 CMD re- 
gions, which efficiently cover the MS, the subgiant branch 
(SGB) and the RGB, and which are large enough to con- 
tain at least a few tens of stars in a typical observed CMD. 
These regions will therefore be able to detect the contri- 
butions from different populations, i. e., stars with differ- 
ent ages and chemical compositions. Figure El shows the 
chosen CMD regions, along with a synthetic CMD made 
up of m = 8 partial models (2 partial models in each 
panel) , whose age ranges are indicated. Most partial mod- 
els are of younger ages. This choice is due to the fact that 
CMDs features are more sensitive to younger turn-offs (1 
- 6 Gyrs) than to older ones. We also tried other sets of 
partial models, with m = 7 and m = 5. In appendix A 
we show the results of several control experiments which 
we used to test the method and its sensitivity to m. As a 
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result of these tests, we chose set A (with m = 8) as the 
most capable of recovering an unbiased SFH. 



5. Results 

The CMD modelling method that we adopted followed a 
regular parameter grid in the SFH: we considered all SFHs 
possible, by varying the (i = 1, 8) coefficients from to 
1, in steps of 0.1, and respecting the constraint given by 
eq. ITI). As for the other model parameters, it was impos- 
sible to fully explore parameter space. Therefore, we kept 
the IMF fixed at all times to that from Kroupa (2002) and 
tested 4 values of E(B-V), 3 values of fbin = 0.25, 0.50, 0.75 
and 3 values of (m — M)o (the one listed in Castro et al. 
2001 and the other two by varying it by ± 0.1). In order 
to account for LMC depth effects, individual model stars 
were assigned (m — M)o values Gaussianly scattered by 
(J(m~M)o — 0.02 around the model distance. This repre- 
sents ~ ±500pc distance spread relative to LMC disk mid- 
point. We considered as our best solution for E(B — V), 
fbin, (ni — M)o in each field the one that resulted in an 
absolute xt minimum in the entire parameter grid. Our 
best SFH solution was then the average values among 
those which satisfied the criterion given by eq. EH We 
tested the sensitivity of the SFH solution (as well as its 
{xt)min value) to E(B — V), distance modulus and fbin- 
The results are shown in §5.4. 

Figures [7| to |51 show the SFHs solutions recovered for 
the different LMC fields, all of them expressed in terms 
of SFH/ < SFH >. As mentioned before, we took the 
average SFR / < SFR > values within the criterion given 
by eq. (|10|l and the la deviations from the average as un- 
certainty bars. A common feature in all recovered SFHs is 
the existence of a non-negligible old (t > lOGyrs) popu- 
lation, as had been previously pointed out by Holtzman et 
al. (1997), Geha et al. (1998), among others. Most SFHs 
deviate significantly from uniformity, displaying periods 
of enhanced star formation {SFR/ < SFR > 2) and 

of quiescence (where SFR/ < SFR > ~ 0.5). 



5.1. NGC1805 and NGC1818 

In figure top panel, we show our best SFH solu- 
tion for the field close to NGC1805, obtained with 
E(B-V)=0.03, fb^n = 0.75 and (m - M)o = 18.59. 
The lower panel shows the result for the field close 
to NGC1818. The absolute minimum (xt)min for 
NGC1818 was found for E(B-V)=0.05, = 0.5 and 
(m - M)o = 18.58. Both SFHs depicted in Figure d are 
quite distinct from the uniform case. Relatively recent 
episodes of star formation have occurred in the past ~ 2 
Gyrs. A broad and significant gap in star formation is 
present in NGC1805 at 6- 10 Gyrs ago. In NGC1818, our 
results point to a star formation closer to uniform before 
2 Gyrs ago; the apparent decrease in SFR in the interval 
6 — 10 Gyrs is still present, but is marginally significant. 



n ' ' ' ' \ ' ' ' ' r 

NGC1805 X„in = 3-23 




Fig. 7. The star formation histories (SFHs) recovered 
from NGC1805 (top panel) and NGC1818 (bottom panel) 
fields, using the E(B-V) values that yield the minimum 
xl (E(B-V) = 0.03 and E(B-V) = 0.05 for NGC1805 and 
NGC1818, respectively). The lines show the average of all 
values of ^ which satisfy the criterion given by eq. (|10() . 
The error bars show the la scatter from the average. 
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Fig. 8. The star formation histories (SFHs) recovered 
from NGC1831 (top panel) and NGC1868 (bottom panel) 
fields, using the E(B-V) values that yield the minimum 
xl (E(B-V) = 0.03 and E(B-V) = 0.00 for NGC1831 and 
NGC1868, respectively). The lines show the average of all 
values of ^ which satisfy the criterion given by eq. (|10|) . 
The error bars show the 1 a scatter from this average. 



5.2. NGC1831 and NGC1868 

The best SFH solutions for the fields close to NGC1831 
and NGC1868 are drawn in Figure|Hl The recovered SFHs 
for both fields are more consistent with uniformity than 
the previous two. The quiescent period between 6 — 10 
Gyrs ago is marginally seen in the NGC 1868 solution 
and is totally absent in the NGC 1831 field SFH. In 
the NGC1868 field an increase in star formation in the 
past T ~ 1 ~ 3 Gyrs is also apparent. The best values 
for distance modulus and binary fraction for this field 
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Fig. 9. The star formation histories (SFHs) recovered 
from NGC2209 (top panel) and Hodge 11 (bottom panel) 
fields, using the E(B-V) values that yield the minimum 
xl (E(B-V) = 0.11 and E(B-V) = 0.10 for NGC2209 and 
Hodge 11, respectively). The lines show the average of all 
values of ^ which satisfy the criterion given by eq. H10() . 
The error bars show the 1 a scatter from this average. 



were (m — M)o = 18.45 and fhin = 0.50, respectively. For 
NGC1831, we found (m - M)o = 18.48 and fun = 0.75. 
The E(B — V) values corresponding to {xVlmm- E(B-V) = 
0.03 for NGC1831 and E(B-V)=0.00 for NGC1868. 

5.3. NGC2209 and Hodge 11 

In figure |51 one can find the results for the NGC2209 
and Hodge 11 fields. For the first we found E(B-V)=0.11, 
(m — M)o — 18.39 and fbin — 0.5 based on the absolute 
iXu)min values, whereas for Hodge 11, our minimization 
procedure yielded E(B-V) = 0.10, (m - M)o = 18.34 and 
fbin — 0.25. The best SFH solutions for these two fields are 
similar. They both point to an increase in star formation 
within the interval 1 — 6 Gyrs ago, especially in the case 
of NGC 2209, where a significant peak is seen at r ~ 1 
Gyr. Star formation seems to have subdued significantly 
since then, according to these SFH solutions. 

5.4. Testing the best solution parameters 

Our method of CMD modelling yields non-unique solu- 
tions for the SFH in a given field; different sets of param- 
eter values often result in acceptable solutions. Figure [TUI 
shows the SFH solutions in the field close to NGC1818 
corresponding to different choices of E(B — V) (panel a), 
(m — M)o (panel b) and /f,™ (panel c). The best val- 
ues of E(B-V), fbin, (m — M)o and the resulting absolute 
ixt)m.in are quoted on the upper panel. In each panel, we 
show the sensitivity of the SFH solution to varying one 
grid parameter at a time. The figure shows that the SFHs 
vary modestly when model parameter values close to the 
best solution are considered. Changes in (m — M)o tend to 
be more important when one wants to find a meaningful 
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Fig. 10. Different star formation histories (SFHs) recov- 
ered from the field near NGC1818, resulting from adoption 
of different E(B — V) values (panel (a)), (m — M)o (panel 
(b)), fbin (panel (c)). (a) The different lines correspond 
to different choices of E(B - V): solid - E(B - V) = 0.05; 
dotted - E(B - V) = 0.03; dashed - E(B - V) = 0.04. 
(b) The different lines show solutions with different values 
adopted for (m - M)o: solid - (m - M)o = 18.58; dotted - 
(m - M)o = 18.68; dashed - (m - M)o = 18.48. (c) The 
different lines show solutions with different values adopted 
for fbin- solid - fbin = 0.5; dotted - fbin = 0.75; dashed 
- fbin = 0.25. In all cases, they correspond to the aver- 
age of all values of ^ which satisfy the criterion given by 
eq. (|10|) . The error bars show the 1 cr scatter from this 
average. 

solution for the SFH. Figure shows {xt)min obtained 
for NGC1818 with the same parameters whose SFH solu- 
tions are presented in figure [TUI One can see the behavior 
of ixDmin with respect to E(B — V) (panel a), (m — M)o 
(panel b) and fbin (panel c) . It is remarkable that a global 
ix'i)min was found in this case. The same behaviour is also 
observed in the other fields. 

As described in the previous sections, our choice of best 
solution was the one which yielded the absolute minimum 
in (xl) with respect to E(B — V), (m — M)o and fbin- We 
also tested our CMD modelling approach by comparing 
our best E(B — V) and (m — M)o to other values of the 
same parameter, found independently by other methods. 
Castro et al. (2001) value of distance modulus was based 
on a model for the LMC disk, as described in §4.1. Their 
extinction values were based on isochrone fits to the CMDs 
of their sample of field LMC stars. Santiago et al. (2001) 
also applied isochrone fits, but to the observed CMDs of 
the nearby rich LMC clusters. Since the clusters are 7.3 
arcmin from the fields studied here, we assume that any 
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Fig. 11. {xi)min values recovered from the field near 
NGC1818, resulting from adoption of different E(B - V) 
(panel (a)), (m - M)o (panel (b)), fbiri (panel (c)). 

Table 3. The best E(B — V) are shown in columns a — d 
while (m — M)o values are listed in columns e, f . (a, e) - 
our results; (b, f ) - Castro et al. (2001) results; (c) values 
quoted by Santiago et al. (2001) for LMC clusters, (d) - 
Burstein & Heiles (1982) values. 



Field 




E(B 


-V) 




(m- 


M) 




a 


b 


c 


d 


e 


./ 


NGC1805 


0.03 


0.04 


0.04 


0.03 


18.69 


18.59 


NGC1818 


0.05 


0.04 


0.03 


0.03 


18.58 


18.58 


NGC1831 


0.03 


0.0 


0.015 


0.01 


18.48 


18.58 


NGC1868 


0.0 


0.02 


0.0 


0.01 


18.45 


18.55 


NGC2209 


0.11 


0.07 


0.11 


0.11 


18.39 


18.39 


Hodge 11 


0.10 


0.0 




0.09 


18.34 


18.34 



variation in extinction across this distance is negligible. 
Notice that this assumption may well be wrong in the case 
of NGC1805, where the broad RGB, seen both in the field 
and cluster CMDs, may indicate differential redenning. 
(m — M)o = 18.50 was assumed in these isochrone fits for 
all clusters. Finally, we also use the E(B — V) values from 
the Burstein & Heiles (1982) maps, which are based on 
galaxy counts, for comparison. 

Table 01 shows our best E(B — V) (column a) and 
(m — M)o (column e) values provided by the partial mod- 
els method. In the other columns we show Castro et al. 
(2001) values for E(B - V) (column b) and (m - M)o 
(column f). The alternative E(B —V) value by Santiago 
et al. (2001) is listed in column c. Both works used Padova 



isochrones from Girardi et al. (2000) in their fits. Burstein 
& Heiles (1982) E(B — V) values appear in column d. 

There is an excellent agreement among the different 
determinations of extinction and distance modulus for 
these fields, especially if we consider the different methods 
and assumptions that are involved. The only exception is 
the very low extinction value adopted by Castro et al. 
(2001) for the field close to Hodge 11, which, according to 
those authors, should be considered cautiously as a possi- 
ble result of a wrong distance modulus or mean metallicity 
adopted for the data. 

Finally, we also carried out the same CMD modelling 
analysis with a cruder version of the chemical-enrichment 
law presented here in order to test the sensitivity of SFHs 
to the assumed chemical enrichment. The SFHs show only 
mild changes but the main characteristics remain, includ- 
ing the existence of field-to-field variations in star forma- 
tion. 

6. Summary and conclusions 

In this work we have analyzed 6 deep CMDs of field stars 
in the LMC using HST/WFPC2 data in the V555 and 
filters. Our fields are located from 1.5° to 6° from the 
LMC centre and within ~ 7' of rich LMC clusters. From 
800 to 3500 stars were measured in both V and I bands in 
each field. We combined these data with those obtained 
by Castro et al. 2001, carefully placing both datasets in 
a homogeneous photometric system, which resulted in a 
much larger photometric sample. 

We have studied the behavior of photometric uncer- 
tainties and of sample incompleteness in the CMDs. We 
developed, tested and applied a method to correct for pho- 
tometric incompleteness as a function not only of magni- 
tude V, but also of (V-I) color. 

Based on the final observed CMDs, and with the aid 
of synthetic model CMDs and objective statistical meth- 
ods, we have constrained the LMC star formation history. 
Our synthetic CMDs are capable of incorporating obser- 
vational effects like extinction, distance modulus, photo- 
metric errors and unresolved binarity. For CMD modelling 
we used the Partial Models Method (Gallart et al. 1999). 
We have tested the method through control experiments 
using 3 different sets of partial models. We also examined 
the sensitivity of our results to adopted distance modulus, 
binary fraction and extinction. 

Although the uncertainty bars are considerable, the 
recovered SFHs vary from one field to another. The two 
fields closer to the LMC bar, NGC1805 and NGC1818, 
show less uniform SFHs, with stronger events of star 
formation. Star formation peaks are seen in the past 
^ 2Gyrs. A relatively quiescent phase from ~ 6 to '--^ 10 
Gyrs ago is also observed, with significant star formation 
earlier on. Periods of enhanced star formation at interme- 
diate ages (~ 2 — 6 Gyrs) are observed in the fields on 
the eastern side of the LMC, NGC2209 e Hodge 11. The 
recent star forming activity (t < 1 Gyr) seen in NGC1805 
and NGC1818 is no longer observed. Finally, a more uni- 
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form star formation history results from CMD modelling 
of the fields in the NW part of the LMC disk, NGC1831 
and NGC1868, especially in this first. The only important 
variation from a uniform SFR in NGC1868 is possibly an 
enhanced star formation period from r ~ 1 to r ~ 3 Gyrs. 
In all fields, one can clearly observe the significant contri- 
bution of older populations (r ~ 6 Gyrs) to their CMDs. 

In brief, we can partially relate the recovered SFHs 
to fields position: fields near the LMC bar (NGC1805 
and NGC1818) present more recent star formation (r 
^ 1 — 2Gyrs) , something that is actually confirmed by a 
larger number of upper MS stars on their observed CMDs. 
NGC2209 and Hodge 11, located to the east of the LMC 
center, show little recent star formation, but display an 
excess of intermediate age (r ~ 2 — 6 Gyrs) stars. Finally, 
NGC1831 and NGC1868 are a further step closer to a uni- 
form SFH, although some increase in star formation in the 
last ~ 3 Gyrs is also visible in the NGC 1868 field. 

It is important to recognize that the inferred SFHs 
do not describe the star formation rates that took place 
in situ, for the stars currently found in our HST/WFPC2 
fields have certainly dispersed away from their various for- 
mation sites. Therefore, the periods of enhanced forma- 
tion we find should refiect star formation rates over much 
broader regions of the LMC, perhaps the entire galaxy in 
the case of the older populations. The same applies to the 
intervals of lower than average activity. The observed dif- 
ferences from one field to the other show that the mixing 
of stellar populations within the LMC field is not 100% 
efficient; on the contrary, there exist coeval stars, result- 
ing from the same burst, which move coherently in the 
galaxy's potential, regardless of their ages or original lo- 
cation. A dynamical model, taking into account typical 
velocity dispersions and orbital motions, is then required 
in order to reconstruct the star formation sites, burst time- 
scales and locations, from our resulting SFHs. 

This work gives continuity to a historical process of 
recognition of the importance of old and intermediate-age 
populations in the LMC. Until the beginning of the 80s, 
the LMC was seen as being formed of relatively young 
stars (Butcher 1977, Stryker 1984). In the beginning of 
the 90s, Bertelfi ct al. (1992) and Vallenari et al. (1996) 
suggested a star formation burst at r ~ 4 Gyr and found 
evidence for some spatial variation in the star formation. 
The HST observations revealed a considerable portion of 
intermediate-age and old stars. The recent reconstructed 
SFHs lead us to conclude that stars older than 4 Gyr are 
at least as numerous as younger ones (Holtzman et al. 
1997, Geha et al. 1998). Our results clearly corroborate 
this notion. 

Another important result is that the gap in star for- 
mation in the LMC is not seen in the SFHs derived here; 
evidence for a decrease in the SFR in the interval 6—10 
Gyrs is seen only in the fields closer to the LMC bar. It is 
therefore hard to reconcile the SFH of fields stars with the 
observed age distribution of LMC clusters, which reveals 
a longer and older age gap. 



It is important to remember that we have always used 
the same chemical enrichment model which is consistent 
with metallicities and ages of LMC clusters (Olszewski et 
al. 1996). However, there is evidence, including from the 
present work, that the field stars SFH differs from the 
SFH of the LMC clusters (Holtzman et al. 1999). Thus, it 
would be extremely important to obtain direct and strong 
constraints on the chemical enrichment of LMC field stars 
(Cole ct al. 2000, Smecker-Hane et al. 2002). For this chal- 
lenge one needs the new generation telescopes, such as 
GEMINI and SOAR. 
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tervals covered by each set. The best solution, as described 
in §5, and which uses the quality criterion given in §4.2, 
is shown as open circles. The open triangles correspond to 
the sets {ai} that yield {xt)mun i-e., the minimum (xt) 
value among all {a^} sets. This {xt)min value is given in 
each panel. The figures show a total of 9 input SFHs, 3 
per figure. 

In most cases, the recovered SFHs (either the best or 
minimum solution) are very similar to the input ones, 
which shows that the method is in fact capable of recon- 
structing a SFH from CMD data. However, some discrep- 
ancies significantly beyond the error bars are seen in a 
few points, especially in panels A.2.C, A.2.f, A. 3. a, A.3.d, 
A.3.f. In most these cases, only one point along the SFH 
is not well recovered, as is the case of the underestimate 
in star formation rate (SFR) at r ~ 8 Gyrs in panel A.2.C. 
This latter panel actually shows a situation (for r ~ 0.5 
Gyr) where the best solution is a good description of the 
input model, whereas the minimal solution is not. This 
is also seen in panels A.3.b and A.3.e. We thus conclude 
that the adoption of an average solution over some range 
in (xt) close to the minimum, as in our best solution, 
rather than the minimum itself, results in more stable re- 
sults. Another conclusion is that SFHs dominated by a 
single burst tend to be more easily reproduced, while more 
complex or nearly constant star formation input scenar- 
ios often result in recovered SFHs with misplaced peaks 
(A. 3. a, A.3.C and A.3.f) or slightly distorted SFH shapes 
(as in A.2.C, A.3.b, A.3.e, A.3.h). Based on the general re- 
sults from these control experiments, and considering the 
relative frequency of discrepancies between input and re- 
covered SFRs, we decided to use set A of partial models, 
with m — 8; this set shows about the same number of dis- 
crepant points as set B (m ~ 7) and set C (m — 5), but 
samples better the SFH. 



Appendix A: Control Experiments 

In order to implement the PMM, we need to define the 
number of partial models to use, m, and to check the effi- 
ciency of the method. This is done by generating an arti- 
ficial CMD based on some known input SFH, as described 
in Sect. 4.1, and submitting this CMD to the modelling 
process described in Sect. 4.2. In other words, we compare 
this "artificial data" CMD with other artificial CMDs cre- 
ated in the same way but following alternative schemes of 
SFH. This comparison makes use of the PMM, with differ- 
ent choices of m. The results of such control experiments 
can be viewed in figures lX!TllA.3l In these figures, the pan- 
els along each horizontal line refer to the results of a given 
set of partial models. We tested three sets: A (m = 8), 
B (m = 7) and C (to = 5). The vertical columns of the 
panels correspond to a given input SFH. The input SFH 
is given as a solid line in each panel. It varies slightly from 
one set to the other because of the variation in the age in- 
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Fig. A.l. Each column shows SFH schemes given in SFR/ < SFR > (solid lines) and corresponding recovered SFH 

(open circles). Each line shows the results for each set of partial models: set A (panels (a), (b) and (c)); set B (panels 
(d), (e) and (f)) and set C (panels (g), (h) and (i)). The set of Oi values obtained on best solution (i. e. {xDmin)- 
Error bars show acceptable a, values with dispersion within la. 
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Fig. A. 2. Same as shown in figure IXm but for the other 3 SFH schemas. 
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Fig. A. 3. Same as shown in figure IXm but for the other 3 SFH schemas. 



